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Partition functions, multiplicity distributions, and isospin fluctuations are calculated for canon- 
ical ensembles in which additive quantum numbers as well as total isospin are strictly conserved. 
When properly accounting for Bose-Einstein symmetrization, the multiplicity distributions of neu- 
tral pions in a pion gas are significantly broader as compared to the non-degenerate case. Inclusion 
of resonances compensates for this broadening effect. Recursion relations are derived which allow 
calculation of exact results with modest computer time. 
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I. INTRODUCTION 

Motivated by the observation of large fluctuations 
of the ratio of neutral to charged particles in cosmic 
ray events H, numerous studies of isospin fluctua- 
tions have been undertaken during the last decade. It 
has been proposed that the melting and subsequent re- 
condensation of the chiral condensate could provide a dy- 
namical means for coherent pion emission where dozens 
of pions are emitted with the same isospin. If N pions 
are confined to a single quantum state in addition to 
being in an isosinglet the probability of finding uq 
neutral pions takes a simple form in the limit of large 
N = n+ + + no, 



dN 



1 



2V7' 



/ = no/N. 



(1) 



One can obtain the same result by considering a coherent 
state 



\ff} = exp(?7 • 7f) |0), 



(2) 



where the pion field operators are ttq = tt^, tt± — 
{ttx ± iTry)/V^, and the direction of ff is averaged over 
all directions. The source of the field ff has been pro- 
posed to be the chiral condensate which might disorient 
itself in a quenching scenario Q. This is often referred 
to as disoriented chiral condensate (DCC). 

The dramatically broad isospin distribution of Eq. (|^) 
relies on the assumption that the emission proceeds via 
a single quantum state. The inclusion of Bose-Einstein 
effects in the thermal emission of pions from a non- 
degenerate array of states was shown to broaden the 
isospin distribution but not nearly as much as in 
Eq. (|l|). Neglecting isospin conservation in accounted 
for the reduced broadening of the peak. A crude account- 
ing for isospin conservation was suggested by considering 
the emission of neutral pion pairs (2/3 tt+tt" and 1/3 
2tto) but came far short of considering the complete 
ensemble of isoscalar states in a multi-level system. The 
emission of pairs through a classical isoscalar field into 



non-degenerate single-particle levels, which may be con- 
sidered as an oriented chiral condensate, has been studied 
as well 0. 

The effects of exact charge conservation in canonical 
ensembles were also studied in other contexts with a pro- 
jection method. Strangeness and baryon number conser- 
vation were found to enchance strangeness productions 
in pN collisions, particularly for small systems The 
confinement of the quark-gluon plasma to color-singlets 
was shown to lead to a reduction in the number of inter- 
nal degrees of freedom, which could lead to measurable 
finite size effects in relativistic heavy ion collisions ||] . 

In this paper, we extend the sophistication of statistical 
treatments by considering the entire ensemble of isoscalar 
states available in a system with many single-particle lev- 
els. Methods for calculating isospin distribution are pre- 
sented, which include Bose-Einstein symmetrization, the 
effects of resonances and the conservation of both total 
isospin / and its projection M. We present sample cal- 
culations to illustrate the effects mentioned above and 
find that symmetrization effects are important for high 
quantum degeneracies, that isospin conservation has lit- 
tle impact when the size of the domain exceeds a dozen 
pions, and that resonances can strongly narrow the dis- 
tribution. 

The next section provides a description of the recur- 
sive methods developed to determine canonical partition 
functions, which include Bose-Einstein symmetrization, 
isospin conservation, and resonances. Techniques for ex- 
act calculation of isospin distributions, which are con- 
sistent with such ensembles, are presented in Sec. Ill, 
results in Sec. IV, and conclusions in ^ 



II. 



RECURSION RELATIONS FOR PARTITION 
FUNCTIONS 

A. Non-degenerate Systems 



When the number of available states is much larger 
than the number of particles, the probability for two 
or more particles occupying the same quantum state is 



2 



small. In such a non-degenerate system, quantum statis- 
tics can be neglected. The partition function for a canon- 
ical ensemble of A particles conserving an additive quan- 
tum number or a vector of such quantities Q can be writ- 
ten as a product of single-particle partition functions 



N 



Za,q = 



(3) 



where N is the number of particle types, Vk is the oc- 
cupation number of particle type fc, is the charge 
of one particle, and the particle number indicates 
how many times a particle contributes to the main con- 
served quantity A. For example, if A is the number of 
pions then ap — 2, since the p meson decays predomi- 
nantly into two pions. The single-particle partition func- 
tion ujk = 5fc exp(— ep'*/T) sums Boltzmann factors 
weighted by the the spin degeneracy gk over all available 
single-particle levels i. 

Summation over the immense number of partitions in 
Eq. (^) can be avoided by rewriting the partition function 
as a recursion relation ^ |2j , 

N 

^A,Q ^^—^ZA-ak,Q-qk- (4) 
k=l 



If intermediate values of the partition function are stored, 
the computations required for Eq. (^) scale linearly in 
both A and TV, thus making it possible to quickly calcu- 
late the canonical partition function numerically. 

A partition function conserving total isospin as well as 
additive quantum numbers is derived by adding a sum 
over all possible isospin configurations for a given parti- 
tion, {vk}, and isospin weights f (/, M\{vk}) to Eq. (0), 



A,I,M 



E 



N 



i = Q) 



(5) 



To convert this partition function into a recursion rela- 
tion insert \Y^k=i '^^i^k = 1 into Eq. (||), 



I,M 



N 

EakWk 
A ^ ivk 



E 



1)! w. 

' J^k J 



J2^{I,M\{,y,}). 



(6) 



The entire system can be broken into two subsystems, a 
single particle with isospin Ik and projection irik and a 
remainder system with isospin /' and projection M — uik, 
which are coupled with the appropriate Clebsch-Gordan 
coefficients. All possible values for the total isospin of 
the remainder system have to be summed over, 



ik+i 

^e(/,Af|{^,}) = 



{Ikmk]I',M-mk\IMf. (7) 



With this modification, the summation indices in Eq. 
can be switched, 



(8) 



and the partition function written as recursion relation 



N ^^^^ I+Ik 

^A,I,M = E~;4~ E ^A^a^.P. 
k=l r = \I-I^\ 



M—mt, 



{Ikmk\l',M ~mk\IMf 



(9) 



Since the partition function is the trace of an isoscalar, 
i.e. e~^l^ , it will not depend on the isospin projection 
M . Hence, Eq. can be further simplified by summing 
the RHS over all isospin projections nik of an isospin 
multiplet. 



^A,7 = ^^^^ E ^A-a^,!', 



I+Ik 



A 



(10) 



where the sum over k' includes iso-multiplets, not indi- 
vidual particles species. 



B. Degenerate Systems 

In a degenerate system several particles might occupy 
the same quantum state, therefore, symmetrization of the 
wave function has to be accounted for. For the purpose 
of this paper we will restrict ourselves to studying Bose- 
Einstein particles. States with multiple particles have to 
be added to Eq. (^) , which only contains states that are 
occupied by zero or one particle. 



oo N 



ji=l k=l 

where the cycle diagram 



-7iak ,M —nrrik i 



(11) 



(12) 



I 



Here, the state \a) refers to an n-particle state of distin- 
guishable particles and \a) is the cyclic permutation of 

(k\ 

that state. The single-particle energy levels are 6; for 
particle type k. A more rigorous derivation of Eq. ( |Tl| ) 
is given in 1 14 . 



Since the partition function for conserved total isospin 
^A,i,M is independent of the isospin projection M, 
as mentioned above, a simple relation can be derived, 
Za,m = X]/>A/ ^A,i, which in turn leads to 



VLa,i = Za,m=i — Za,m=i+i- 



(13) 
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A second method for calculating the pion partition 
function constraining total isospin is obtained by eval- 
uating the cycle diagram Cn'^'* in Eq. (|l2|) for its isospin 
content. These new cycle diagrams are defined as 



C„,. = ^(/3,n,z|e-^/^|/3,n,i 

/3 



(14) 



where the sum over /3 represents a sum over all states with 
fixed particle number n and isospin i. The particles are 
assumed to be distinguishable and P represents a cyclic 
permutation of particles. The partition function for the 
pions in term of this new cycle diagram is then 

^ A n I+i 
^^^■^"'tX^X! X! Cn,i^^A-n,/', (15) 



=1 i=o r=\i-i\ 



where the new cycle diagrams C are yet to be determined. 
After obtaining a recursion relation for these functions 
from Eq. (Ill) itself, 



A-l n I+i 

CAJ^AnAJ-J2Yl Cn,^^A-nJ', (16) 

n=l 1=0 I' = \I-i\ 

we find that these cycle diagrams follow a simple pattern 
by considering a one-level system where il,Aj is easily 
calculated. 



Cn,i 



Cm i — n 

Cm i^O 
0, otherwise 



(17) 



Since resonances are more massive and have lower 
phase space occupations, the probability of creating sev- 
eral resonances in the same state can be neglected except 
in the limit of extremely high densities. Therefore, reso- 
nances might be treated as independent, non-degenerate 
subsystems, for which a partition function can be ob- 
tained through Eq. (^ . The partition functions of two 
subsystems, 1 and 2, can then be convoluted to obtain 
that of the entire system. 



A,I,M 



A A' 

EE 

A'=OI'=OI' 



I+I' 

E 



A', I" 



(18) 



If there are more than two subsystems Eq. (|18|) can be 
applied successively, i.e., partition functions of any two 
subsystems are convoluted first to obtain a new parti- 
tion function, which is then convoluted with the partition 
function of another subsystem, and so on. 

It should be pointed out that treating resonances as if 
they are in a different system is not consistent with the 
indistinguishability of pions from resonances and direct 
pions. For narrow resonances, one could uniquely iden- 
tify the pions by constructing the invariant masses of the 
constituents. However, for broad resonances, one can not 



confidently identify the resonances. The criteria for res- 
onances being narrow are identical to the criteria that 
they have sufficiently long lifetime to decay outside the 
collision region. For all the calculations considered in this 
paper, it is assumed that the resonances are separable. 
This assumption is excellent for pions from rj mesons, 
good for pions from uj resonances and questionable for 
pions from p decays. Including the effects of symmetriz- 
ing the resonant and non-resonant pions remains an open 
question. 



III. MULTIPLICITY DISTRIBUTIONS AND 
ISOSPIN FLUCTUATIONS 

A. Non-degenerate Systems 

The multiplicity distribution can be calculated from a 
ratio of partition functions, where the numerator includes 
an extra constraint, 



PA^Aiirij) 



Za. 



Ma 



Za. 



M 



(19) 



Here the numerator represents a canonical ensemble with 
the appropriate conservation laws containing rij particles 
of type j. This additional constraint can be regarded as 
a "charge" and, therefore, is added to the indices of the 
partition function. The partition function in the numer- 
ator can be rewritten with the aid of Eq. (H), 



Za. 



Ma 



N 
k=l 



A 



ZA-au 



(20) 



where the feed-down factor dk.j indicates that a parti- 
cle of type k decays into dkj particles of type j. This 
approach, however, will not work for a resonance that 
decays via more than one decay channel. In such a case, 
a pseudo-particle is included for each decay branch with 
its degeneracy gi scaled by the corresponding branching 
ratio. 

It will prove more convenient to write equations in 
terms of the product of partition function and multiplic- 
ity distribution 



WA,M{nj) = ZAAiPAAiinj) 



(21) 



instead of the multiplicity distribution itself. For a non- 
degenerate system conserving only additive charges, the 
multiplicity distribution can be obtained from 

N 

WA,Minj) = -^^A^a^M-mA^j - dk,j). (22) 
fc=l 

The occupation number for particle type j including 
feed-downs from resonance decays is determined by mul- 
tiplying the occupation number of all particles, given in 



4 



[ [12[ , by the feed-down- factor and summing over all reso- 
nances, 



1 ^ 



(23) 



The second moment of the distribution will be needed for 



calculating isospin fluctuations as in Sec. [II C. 

Edk.jdk' .j' r e „ 
\Ok,k'UJk^A-ak,M-m.k 



M 



k,k' = l 

}. (24) 

The multiplicity distribution incorporating conserved 
total isospin can be derived in a similar manner as was 
employed for Eq. 



N 



W^A,/,M(n,) = ^-^^(/fcmfc;/',Af-mfc|/Af)2 
fc=i /' 



W^A-a,,/',M~m,('^j - dk,j). (25) 



B. Degenerate Systems 

When considering only quantum numbers correspond- 
ing to additive charges, the multiplicity distribution is 
obtained through Eqs. (|l]) and ([l9| ) and we find 

N oo 

Wa,m{'^3) = X! ^ X! ^I'^^^A-lajM-lTyi, [rij ~ ldk,j), 
k=l 1=1 



(26) 



where Ci is the cycle diagram define in Eq. (|12|). Cal- 
culating the two-point function and the four-point func- 
tions allows one to obtain the first two moments of the 
multiplicity distribution. The 2-point function is 



Za. 



M 



exp(-nei/T)ZA-„Q, 



M- 



(27) 



Summing over all particle types and states and multiply- 
ing by feed-down factors results in an expression for the 
occupation numbers, 



N oo 



k=i 1=1 
Similarly, the 4-point function 

ft „ \ _ ^U^jk + Sik^jl 



k ,M~^'^k ■ 



(28) 



{a]a]akai) 



Za,m 

exp{-niejT) ex.p{-njej /T) 



ZA—Uiai 



(29) 



serves to obtain second moments of the multiplicity dis- 
tribution 



Z 



^'^ Cij^i'Za- 



A.M 



LI' 



Za. 



I y^,Sj.j'(^lZA-la 



y^ ClCl' ZA-lak-l'a^, 



k,M — lmk 



Al — lrUk —I' ray 



LI' 



(30) 



Neglecting the terms with Cg where i > 1, one returns 
to the non-degenerate result, Eq. (|2^). 

Calculating the multiplicity distribution for degenerate 
systems with the constraint of total isospin conservation 
becomes difficult because the analog of the cycle diagram, 
needs to be analyzed for isospin i and 
1 a cycle diagram with a particles and 



Ci^^ in Eq, 



charge nj. Sue 
isospin projection m can be written as 

Xa,i,rn{nj) = ^{ct, a, i, m\e^ " ^'^ \/3 , a, Tij , m) (31) 

•(/3, a, rij, m\a, a, i, to), 

where the sums over a and f3 correspond to sums over 
all states with fixed {a,i,m) and {a,nj,m), respectively, 
and a represents a cyclic permutation of particles, which 
are assumed to be distinguishable. 

The multiplicity distribution can be calculated in 
terms of these cycle diagrams 



WA,I^M{nj) = ^ ^ Xa,i,m{n'j) 



(32) 



a.i.m.n ■ 



■WA~aj',M-minj ~ n'^)(l', M - to; i, m\I, 

To derive the cycle diagrams x, consider a single-level 
system, for which the probability distribution W^^^ is 
calculated in Appendix It follows from Eq. ( ^^ that 

a<A.i.'m,nj ,/' 

(33) 

The function x which accounts for all levels can be easily 
generated from x'"^\ 

Eff^e-^'/^, (34) 

e 

where £ indicates the single-particle energy levels with 
energy Eg. 

If different particles, like resonances, are to be included 
in the ensemble, one can either insert a sum over species 



5 



into Eq. (^2|); or calculate W separately for each species 
and convolute them to find W for the entire system, 



(35) 



A'I'.M' 



w 



T-A-,i",M-M'inj ~ n'^Hl'M'; I", M ~ M'\IM)' 



In our calculations, W was calculated separately for res- 
onances neglecting symmetrization and then convoluted 
with W calculated for pions with proper symmetrization. 



C. Isospin Fluctuations 

Given some system with quantum states a, pion 
isospin fluctuations can be defined as 



^^{a\{N+ + N_ - 2iVo)^|a) 



(36) 



where Nj^ , N_ and Nq are the number operators of the 
respective pions. These isospin fluctuations could be 
computed through multiplicity distributions or with the 
expressions for densities and higher moments that were 
given above. However, when total isospin should be con- 
served as well, multiplicity distribution calculations are 
slow and expressions for densities and higher moments 
are difficult to derive. Instead, we will derive the isospin 
fluctuation for a system in an isosinglet in terms of isospin 
projection states. 

The operator in Eq. (^ is a product of two rank-2 
spherical tensors components 



{N+ +N^-2Nof = 6T20T, 



20, 



where 

T2±l 
T2±2 



(37) 

(38a) 

(38b) 
(38c) 



A product of spherical tensors can be decomposed into 
other spherical tensor components 



T2oT2o = ^(20;20|JAf)A 



JM- 



(39) 



By the Wigner-Eckart theorem only Aqq contributes 
when contracted between isosinglet states, and because 
^00 is an isoscalar we can write 



6 

V5 



^(Aoo) = 4 ( E - E ) (^-^l^oolM). 



(40) 



After some algebra we obtain 

^00 = E<2Af; 2, -M\m)T2MT2.-M 



A' 



00 



AT, 



M 



(41) 



where 



A' 



1 

V5 



iV_ - 2iVo)' 



(42) 



and 



A, 



QM 
00 



(43) 



The expectation of is non-zero when particles of 

different charges are in the same quantum state, or when 
two differently charged pions are produced into two dif- 
ferent states with a quantum correlation arising from a 
resonance decay. The contribution to from the de- 

generate nature of the pion states can be determined via 
Eq. (H), 



Za. 



M 



Y^illZA- 
l.V 



lai^—l'a^i .M — Imi^—V m^t 



(44) 



This contribution can be ignored in the non-degenerate 
limit, where occupation numbers are small. 

Contributions to from the coherent correlation 

between pions from resonant decays can be found by ex- 
pressing the resonances in terms of pion creation opera- 
tors. For example, the meson, which is a member of 
an isotriplet, can be considered as one pion in an s wave 
and a second pion in a p state. Referring to these two 
states as i and j, 

1 



IP 



x/2 



|0). 



(45) 



How the states i and j are chosen is irrelevant since they 
are summed over in Eq. (^), but the coherent mixture 
of the two permutations, which is necessary for the p"*" 
to be a member of an isotriplet, results in a non-zero 
contribution to ASJ^, 



^00 ' 



{AT) 



1 



V5 



Np+ 



(46) 



The LJ and 77 mesons are isosinglets and can be treated 
accordingly. For instance, 



which adds another term to ^[j^o'^, 



{aT) = -^^{n. 



(47) 



(48) 
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IV. RESULTS 

In heavy ion reactions, and perhaps in pp reactions, 
pions reinteract with other pions in their neighborhood, 
or domain, and might be expected to sample a large por- 
tion of the the available phase space. As isospin should be 
conserved in each domain, it seems reasonable to explore 
distributions for a few dozen pions rather than creating 
an ensemble of a few thousand pions, which could be 
treated as a grand canonical ensemble |13 . One of our 



study's goals is to understand how many pions are re- 
quired for conservation constraints to become irrelevant. 

In the following, symmetrization and resonances are 
first ignored in order to focus on the effects of conserving 
isospin, subsequently, the effects of symmetrization and 
resonances are illustrated by considering a simple exam- 
ple. 



A. Total Isospin Conservation 

When quantum degeneracy and resonances are ig- 
nored, isospin distributions are independent of energy 
levels or temperature. Therefore, the results presented in 
the following are generic to any system where only pions 
are considered and the phase space occupation numbers 
are small. A random distribution ignoring isospin con- 
servation, i.e. a mixed-event construction, will serve as a 
benchmark. 



^random('T'o) — ( g 



N 



E 



i^-\-n_-\-7iQ—N 



n+!n_!no' 



(49) 



Unlike distributions that conserve isospin, this distribu- 
tion allows both even and odd numbers of neutral pions 
and is, therefore, scaled by a factor of two to compare 
the width with that of the other distributions. Secondly, 
when pion creation is constrained to isoscalar pairs, as in 
1^ , the distribution can be considered as a binomial dis- 
tribution of pairs where one third of the time the pair is 
comprised of two neutral pions and two thirds of the time 
the pair is comprised of a positive and negative pion. 



pairwisc(^o) — [ ^ 



no/2 



I) 



(Af-no)/2 



iN/2)\ 



(50) 



(no/2)!(Ar/2-7io/2)! 

This pairwise distribution is broader than the random 
distribution by a factor of \/2, as can be seen in Fig. |^. 

Finally, the isospin distribution for non- degen erate 
particles is calculated with the methods of Sec. [II A| with 
all 12-pion isosinglet states being considered. The con- 
straint of exact isospin conservation only modestly broad- 
ens the distribution relative to the random distribution 
as shown in Fig. 

These findings are underscored by comparing the 
isospin fluctuations as a function of total pion number. 
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FIG. 1: The probability of observing no neutral pions is 
shown for a system of 12 pions. Symmetrization and reso- 
nances are neglected. A random distribution, multiplied by 
2 for comparison, is represented by a solid line and is used 
as a benchmark. A distribution resulting from isoscalar pairs 
of pions (squares) is significantly broader, whereas a distribu- 
tion including all isosinglets (circles) approaches the random 
distribution. 
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FIG. 2: The fluctuation for non-degenerate systems as 
a function of system size for an ensemble of isoscalar pairs 
(squares) and one including all isosinglet states (circles). The 
fluctuations have been normalized by the fluctuation for a 
random system. The fluctuations for the pairwise distribution 
are twice the fluctuations for a random distribution, whereas 
including all isosinglet states relaxes the constraint and re- 
sults in the same width as the random distributions for large 
systems. 



as shown in Fig. ^. When pion emission is constrained 
to isoscalar pairs, fluctuations are twice as large as com- 
pared to the random case for all system sizes. When 
considering all A^-pion isosinglets fluctuations are larger 
by a factor which falls from two to unity as A^ approaches 
infinity. 
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B. Including Symmetrization and Resonances 

To illustrate the effects of Bose-Einstein symmetriza- 
tion, total isospin conservation, and resonance decays, 
an assumption must be made about the available single- 
particle energy levels that are summed over in the single- 
particle partition function. Assuming the model system 
is confined to a cube of volume V, the energy states are 
obtained 



i?2 



(51) 



where R = V^^"^, M is the mass of the particle, and 
n, m, I are chosen to be half- integers. The choice of half 
integers, instead of the more usual integers, deemphasizes 
zero-point surface energy effects and seems more physical 
if the confinement to the volume does not arise from an 
infinite potential well. This becomes important when 
systems are confined to a small volume. 

Isospin distributions were calculated for a system of 
24 pious at two densities, 0.3 fm~^, which is well above 
breakup densities for hadronic collisions and 0.1 fm~'^. 
The temperature was chosen to be 125 MeV. The Bose- 
Einstein nature of pions was taken into account when 
calculating multiplicities according to the formalism from 
The multiplicity distributions were then con- 



Sec. [II B 



voluted with those of three resonant states, the isotriplet 
p mesons and the isosinglet u and rj mesons, which 
were treated as non-degenerate systems according to 



Sec. [II A. In principle, strange mesons and baryonic res- 



onances can be incorporated as well, but the calculation 
would be significantly lengthened by the inclusion of ex- 
tra indices. Prospects for such calculations are discussed 
in the conclusions. 

As expected, the isospin distributions for symmetrized 
pions in an isosinglet are broader than the random distri- 
bution when resonances are neglected, as shown in Fig. ||. 
This broadening is especially pronounced at high density. 
However, the inclusion of resonances more than compen- 
sates for the symmetrization effects and results in distri- 
butions that are narrower than the random distribution. 
Figure ^ displays fluctuations as a function of density for 
the 24-pion system. The dramatic broadening induced 
by symmetrization at high density is counteracted by a 
remarkable narrowing when resonant states are consid- 
ered. 

The conditions for pions to prefer forming a resonance 
are similar to the conditions for symmetrization to be 
important, i.e., a high phase space density. At low tem- 
perature, where the mass penalty for resonant formation 
would play a larger role, resonant effects would become 
relatively less important than symmetrization. However, 
for the temperature of 125 MeV considered here, the res- 
onant effects overwhelm the effects of symmetrization at 
all densities. 

Finally, it should be noted that the isospin fluctuation 
6*2 was calculated both from the distributions themselves 
and from the methods described in Sec. IIIC. Although 
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FIG. 3: The probability for producing no neutral pions in 
a system of A = 24 pions at T = 125 MeV is shown for 
ensembles restricted to an overall isosinglet and with sym- 
metrization included. Calculations with resonances (squares) 
and without (circles) are displayed. Although symmetrization 
broadens the distribution relative to the random distribution 
(solid line), the inclusion of resonances results in a narrower 
distribution as compared to the random distribution. 



the two sets of numerical calculations have little in com- 
mon aside from the functions used to generate the single- 
particle levels, the two sets of moments agreed with each 
other within the numerical accuracy of the computer. 



V. CONCLUSIONS 

For the first time, expressions were obtained for multi- 
plicity distributions and isospin fluctuations for a canon- 
ical ensemble, in which total isospin as well as additive 
quantum numbers are exactly conserved. The formalism 
has been extended to include both pions and resonances 
and can account for Bose-Einstein symmetrization of the 
pion wavefunction. Numerical calculations were then 
performed to study the effects of total isospin conserva- 
tion, quantum symmetrization, and resonance decays on 
the width of the multiplicity distribution, which can be 
squared to obtain the isospin fluctuations. Direct expres- 
sions for the isospin fluctuations permit a much quicker 
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FIG. 4: Isospin fluctuations, or squared width of the isospin 
distributions, as a function of density scaled by the width of 
a random distribution. Calculations were performed for 24 
pions restricted to isosinglet states at a temperature of 125 
MeV. At high density, the distributions are broadened by in- 
cluding symmetrization and narrowed by including resonances 
produces. 



calculation than those for the multiplicity distributions. 

It was found that conservation of total isospin and its 
projection has little effect on the width of the multiplicity 
distributions, when the systems are larger than a dozen 
particles. It should be noted that conserving only the 
projection, not the total isospin, would result in distri- 
butions where the average number of neutral pions would 
not equal one third of the total. At high phase space den- 
sities, including Bose-Einstein symmetrization leads to a 
multiplicity distribution that is much broader than a ran- 
dom distribution. However, addition of resonances more 
than compensates for this broadening and narrows the 
multiplicity distribution below the width of the random 
distribution. Both effects are small when the phase space 
density is below 0.1 fm~^. 

The widths of the multiplicity distributions are largely 
dominated by the behavior of the tails, thus making it 
imperative to perform exact calculations. These calcu- 
lations were made possible by using recursion relations 
that circumvent summing over the immense number of 
partitions in the partition functions. 

The calculations in this paper were largely schematic 
and included only three meson resonance, p, 77, and 
uj. For a more realistic calculation more resonances like 
strange mesons and baryonic resonances need to be added 
along with strangeness and baryon number conservation. 
The formalism presented in this paper scales linearly with 
the number of particle species if the resonances were to 
be included without additional conserved charges. The- 
oretically, any amount of quantum numbers can be con- 
served as long as they commute with the isospin op- 
erator. Practically, every new index, which has to be 
added to the partition function and multiplicity distri- 



butions, to conserve another charge increments the num- 
ber of loops by one. The increase in runtime has little 
consequence for calculations of partition functions and 
direct calculations of isospin fluctuations, which are vir- 
tually instantaneous, whereas any additional index would 
significantly slow down calculations of multiplicity distri- 
butions, which take on the order of ten minutes. 

The aforementioned caveats are not expected to be- 
come major obstacles in the application of the presented 
formalism. Although the particle multiplicities are high, 
in the thousands, in possible physical applications like 
relativistic heavy ion collisions, the local nature of charge 
conservation would limit the number of particles consid- 
ered at any given time. The system under considera- 
tion would have to be broken into domains, in which the 
charges are conserved locally, then calculations proceed 
one domain at a time. Each domain would have a rela- 
tively small number of particles with which to cope. 

Before tackling the more numerically challenging prob- 
lem of including strangeness and baryon number, one 
should consider the limitations of any comparison with 
experiment. Most importantly, it difficult to count neu- 
tral pions as each neutral pion decays into two photons. 
Furthermore, one should consider the ability to identify 
neutrons and kaons, especially those kaons which then 
decay into pions. Given the inherent complexity of any 
such measurement, we felt that it was proper to stop 
short of performing more complicated calculations with- 
out a commensurate consideration of the measurement. 
Nonetheless, several valuable lessons were gained from 
the schematic calculations presented here. 

Another possible application of partition functions 
with exact quantum number and isospin conservation, 
one that has not been explored in this paper, are Monte 
Carlo algorithms for particle generation. Recently, there 
has been much interest in modeling relativistic heavy ion 
collisions with hybrid models, in which early, dense stages 
of the collision are described by a hydrodynamical model 
before switching to a hadronic cascade to simulate the 
freeze-out stage The change of degrees of freedom 

at the interface between the two models from the energy- 
momentum tensor to hadrons is generally modeled by a 
grand-canonical ensemble, which conserves charges only 
in the average over many events. However, event-by- 
event charge conservation is essential to calculating ob- 
servables like fluctuations and balance functions, which 
have been proposed as possible signal for the quark- 
gluon-plasma [pq. 
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APPENDIX A: SINGLE-LEVEL MULTIPLICITY 
DISTRIBUTION CONSERVING TOTAL ISOSPIN 

By definition, the distribution W^Vm('^j)> needed in 
Eq. (^), can be written as the product of the partition 
function and the multiphcity distribution, 



W^^^ =uj'-'-'- v^'' in) 



(1) Jl) 



(Al) 



where the partition function for a particles in a single 
level with energy E is 



eyi^[—aE /T) , if a + i even 
"a,i,ni - I , if a + i odd 



,(1) 



(A2) 



Pal\ m ("j ) is probability of observing rij pions in a 
single-state system that contains a total of a pions with 
total isospin i and projection m. This probability dis- 
tribution has to be calculated for only one type of pions 
because the pion occupation numbers are related through 



a = no + n_|_ + n_ 



(A3) 
(A4) 



In the following, probability distributions will be derived 
for positive pions. 

The isospin wave function of the system can be written 



in terms of eigenstates of the number operators 

(a+m)/2 

|no)|n+)|n_), (A5) 

n-i- — m 

where Hq — a + m — 2n+ and n_ = n+ — m. The coef- 
ficients in Eq. (|A^ ) are related to the probability distri- 
bution by 



(A6) 



The isospin wave function |a,a, a) can only be con- 
structed if all pions in the state are positive, i.e., n_|_ = a, 
therefore 



1 , if n+ = a 
, otherwise 



(A7) 



Leaving the pion number a and isospin i ~ a fixed, we 
can apply the isospin lowering operator /_ = \/2(ji'^_'Kq + 



Tr^TT^) to reach lower values of m, 



I-\a,i,m) = + 1) — m{m — l)\a,i,m — 1). (A8) 
The LHS of Eq. (^) expands to 



I-^aa,a,m.n+\a + m - 2n+)\n+)\n+ - m) 

(a+m)/2 

= %/2 ^ y^a + m- 2n+y^n+ - m + laa,a,m,n+|a + (m - 1) - 2n+)|n+)|n+ - (m - 1)) 

n-l- — m 
(a+T?i)/2-l 

+ \/2 ^ \/a + (m - 1) - 2n+^n+ + laa.a,m,n++i\a + (m - 1) - 2n+)|n+)|n+ - (m - 1)). (A9) 

n-i-— m — 1 



When the coefficients on the LHS are match with those on the RHS of Eq. (A8) a recursion relation is obtained, 



^a,a,m— l,n-|- 



i{i + 1) — m{m — 



Yy I V(a + m - 2n+)(n+ - m + l)aa,a,m,n+ 



+ (a + m - 1 - 2n+)(n+ -I- l)aa,a,m,n+ -i-i f- (AlO) 

I 



So far, we have only found the coefficients aa,i,m,n+ for where the normalization constant is 
a = i. With the help of the isoscalar operator 

U2 = 27r^^7rl - nlni (All) 

that creates two pions without altering the isospin, higher 
values of a can be reached. 



\/(a + 2)(a + 3)-i(i+l). (A13) 



U2\a,i,m) = Na,i\a + 2,?,m), 



, , Matching the coefficients on both sides of Eq. ( A12 ) to 
^ ' each other leads to 
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aa+2 ' 



n' 



>/(a + m — 1 — 2n+)(a + m — 2 — 2n^)aa.i.i 



(A14) 



In case, a + z is odd, no combination of pions yields 
the isospin wave function \a,i,m), as is evident from 
Eq. (^). Therefore 



Equations (jA^), ( |A10| ), ( [Al^ ), and ( [AT^ ) completely 



determine the coefficients 



.(1) 



, which in turn define 



a, 



a,z,m,n4 







if a + z odd. 



the probability distributions Pa imi''^+) ^'^d, therefore, 

.(1) 



(A15) the single-level partition functions 
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